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We present general relativistic correction terms appearing in Newton's gravity to the second-order 
perturbations of cosmological fluids. In our previous work we have shown that to the second-order 
perturbations, the density and velocity perturbation equations of general relativistic zero-pressure, 
irrotational, single-component fluid in a flat background coincide exactly with the ones known in 
Newton's theory. We also have shown the effect of gravitational waves to the second-order, and pure 
general relativistic correction terms appearing in the third-order perturbations. Here, we present 
results of second-order perturbations relaxing all the assumptions made in our previous works. We 
derive the general relativistic correction terms arising due to (i) pressure, (ii) multi-component, 
(iii) background curvature, and (iv) rotation. In case of multi-component zero-pressure, irrotational 
fluids under the flat background, we effectively do not have relativistic correction terms, thus the 
relativistic result again coincides with the Newtonian ones. In the other three cases we generally 
have pure general relativistic correction terms. In case of pressure, the relativistic corrections appear 
even in the level of background and linear perturbation equations. In the presence of background 
O |. curvature, or rotation, pure relativistic correction terms directly appear in the Newtonian equations 

of motion of density and velocity perturbations to the second order. In the small-scale limit (far 
inside the horizon), relativistic equations including the rotation coincide with the ones in Newton's 
gravity. All equations in this work include the cosmological constant in the background world 
model. We also present the case of multiple minimally coupled scalar fields, and properly derive 
the large-scale conservation properties of curvature perturbation variable in various temporal gauge 
conditions to the second order. 
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I. INTRODUCTION 

Large amount of cosmological data on the large-scale structures and motions of galaxies [3, 0] , and the temperature 
and polarization anisotropies of cosmic microwave background radiation 0, Q| have been accumulating recently. In 
current standard cosmological scenario such structures are explained as small (linear) or large (nonlinear) deviations 
from spatially homogeneous and isotropic Friedmann background world model. In order to explain these data theo- 
retically, researchers rely on the linear perturbation theory based on relativistic gravity and quasi-linear perturbation 
theories and nonlinear simulations based on Newton's gravity. To the linear order in perturbation the general rel- 
ativistic result was first derived by Lifshitz in 1946 and later shown to coincide with the Newtonian result in a 
zero-pressure medium The same is also known to be true for the background world model. That is, the general 
relativistic result was first derived by Friedmann in 1922 and later shown to coincide with Newtonian result in a 
zero-pressure medium Q. 

The observed large-scale distribution of galaxies shows that in the largest observed scale (say larger than several 
hundred mega-parsec scale) the distribution may not be inconsistent with the linear assumption around the Friedmann 
background. However, as the scale becomes smaller the distribution apparently shows quasi-linear to fully nonlinear 
structures. The fully nonlinear processes occur in small scale where the relativistic effects characterized by GM/rc 2 ~ 
v 2 /c 2 are quite small. If we could ignore such relativistic effects, Newton's gravity would be sufficient to handle 
the relevant nonlinear processes. If we need to consider the weakly relativistic correction terms in fully nonlinear 
stage, instead of the relativistic perturbation approach which can handle the fully relativistic processes under weakly 
nonlinear assumption, we can use the post-Newtonian approximation developed in the context of cosmology in Q. 

For structures in the quasi-linear evolution phase, previous researches were based on Newton's gravity especially 
assuming the single component zero-pressure fluid without rotational perturbation In our previous works in 
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[lll.ll2j| we have shown that, in the single component zero-pressure fluid without rotational perturbation, cosmological 
scalar-type perturbation equations in a spatially flat background coincide exactly with the Newtonian ones up to the 
second order in perturbation. In [111, E3 we a l so have shown the contribution of gravitational wave perturbations to 
the hydrodynamic parts in the second-order perturbations. In Newton's gravity the hydrodynamic equations of zero- 
pressure fluid contain only the quadratic order nonlinearity. In [l3| we presented pure general relativistic correction 
terms appearing in the third-order perturbation, and showed that all third-order correction terms are 1CP 5 times 
smaller than the second-order relativistic/Newtonian terms, and independent of the horizon scale. 

In this work we will take into account of the pure general relativistic effects appearing in the second-order perturba- 
tions which were ignored in our previous work in [12| ■ We will consider general relativistic effects of (i) pressure, (ii) 
multi-component, (hi) background curvature, (iv) rotation in cosmological fluids to the second-order perturbations. 
As results we will show that, although in [l2j we have shown the exact relativistic/Newtonian correspondence to the 
second-order perturbations by ignoring the above four conditions and the gravitational waves, as we take these four 
effects into account we often encounter pure general relativistic effects appearing in the corresponding Newtonian equa- 
tions even to the second order in perturbations. Our results will show that the relativistic/Newtonian correspondence 
continues even in the multi-component situation assuming zero-pressure irrotational fluid in a flat background, but 
in the presence of the cosmological constant. This is a practically useful result because the matter content of present 
universe is dominated by collisionless dark matter and baryon both of which practically have zero-pressure. In the 
other three cases, relaxing any of the assumptions about pressure, rotation, and background curvature generally leads 
to pure general relativistic correction terms to the second order. We will present such correction terms in the context 
of Newtonian hydrodynamics. One additional relativistic/Newtonian correspondence occurs in the case of rotation in 
small-scale (sub-horizon-scale) limit which is another practically important result. This correspondence allows us to 
use the Newtonian equations safely in such a small-scale limit even in the presence of rotational perturbation to the 
second order. 

In Sec. |TT] we summarize Newtonian hydrodynamic perturbation equations valid to fully nonlinear order. In Sec. 
Mil we briefly summarize our previous result of relativistic/Newtonian correspondence to the second order, and pure 
general relativistic correction terms appearing in the third order. In Sec. IIVI we present parts of the covariant and 
the ADM (Arnowitt-Deser-Misner) equations which are valid in multi-component situation. In Sec. [V] we present 
the basic perturbation equations valid to second order. In (llj the basic set of equations was presented using fluid 
quantities based on the normal-frame four-vector. The fluid quantities in the present work are based on the energy- 
frame four-vector, and in this section we present the basic equations using such fluid quantities. In Sees. IVHIIXI we 
analyse the effects of the pressure, the multi-components, the background curvature, and the rotational perturbation, 
respectively. In Sec. Ewe- present equations in the scalar fields and generalized gravity theories using the energy-frame 
fluid quantities. In Sec. IXII we properly derive conservation properties of curvature perturbation in various temporal 
gauge (hypersurface) conditions to the second order in perturbations. Section IXIII is a discussion. In this work we 
follow notations used in [ll|, [l^] ■ We set c = 1 , but when we compare with Newtonian case we often recover the speed 
of light c. 



II. NEWTONIAN NONLINEAR PERTURBATIONS 



In order to compare properly the relativistic results with the Newtonian ones, in this section we summarize the 
Newtonian cosmological perturbation theory in fully nonlinear context. We consider multi-component fluids in the 
presence of isotropic pressure. In case of n-fluids with the mass densities ft, the pressures pi, the velocities 
(i = 1, 2, ... n), and the gravitational potential <i>, we have 

ft + V-(ft Vl ) = 0, (1) 

v« + Vi • Vvi = Wpi - V$, (2) 

n 

V 2 $ = 4^G]Tft, (3) 

j=i 

Assuming the presence of spatially homogeneous and isotropic but temporally dynamic background, we introduce 
fully nonlinear perturbations as 

Qi^Qi + Sgi, p.^p. + Spi, v l= Hr + Ui, $ = $ + <5$, (4) 

where H = a/a, and a(t) is a cosmic scale factor. We move to the comoving coordinate x where 

r = a(t)x, (5) 
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thus 



V = V r = -V x , 

a 

d__d_ 
dt~ dt 
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X dt 



Hx ■ V x . 



In the following we neglect the subindex x. To the background order we have 



Qi + 3Hg t = 0, 



8itG ^ 2E 



(6) 



(7) 



where E is an integration constant which can be interpreted as the specific total energy in Newton's gravity; in 
Einstein's gravity we have 2E = —Kc 2 where K can be normalized to be the sign of spatial curvature. To the 
perturbed order we have [bj ] 



-V-U; = -~V- (^tO, 

a a 



a 



1 v * w 



agi 1 + Si a 



(8) 
(9) 
(10) 



By introducing the expansion 6i and the rotation a?; of each component as 



1 -> 1 

-V-Ui, W,E-VxUj 

a a 



Eq. ([9]) gives 



~UJi + 2H~LJi = -V X (u,; • VUi 

a 1 



1 (V*<) x VS Pl 



\1 + Si 



a 2 Qi (1 + Si) 2 

By introducing decomposition of perturbed velocity into the potential- and transverse parts as 



Ui^Vui + u^, V.uf»E0; 9 i = ^m, tf, = I V x uj" \ 



instead of Eq. (fTJ]) we have 



u) ' + Hu) ' 



Ui ■ Vu t + - VA _1 V • ( Ui • Vu, 

ft 1 + Oj 



1 VS Pi 
Qi 1 + 



(11) 

(12) 
(13) 

(14) 

(15) 



We note that the pure u% contributions in the right-hand-side of Eq. (| 13[) or Eq. (flSI) vanish. Thus, under van- 
ishing pressure, pure irrotational perturbation cannot generate the rotational perturbation. Equation (| 1 3[) shows 
that presence of pressure perturbation oblique (i.e., non-parallel) to the density perturbation can generate rotational 
perturbation. Combining Eqs. . (fT!?|) we can derive 



5i + 2HS t - 4irGj2 SjSj = ~ [aV • (5iUt)\ + -^V • (iij • Vu;) + -^Y 



a 2 Qi 



1 + 5* 



(16) 



Equations (|8|l- (jT6|) are valid to fully nonlinear order. Notice that for vanishing pressure these equations have only 
quadratic order nonlincarity in perturbations. 
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III. SUMMARY OF PREVIOUS WORK 



In [111 [l2| we have derived second-order perturbation equations valid for the single component, irrotational, and 
zero-pressure medium in zero-curvature background. These are 



5+ -V-u = --V • (Su) , 
a a 



-V • (u + Hu) + AkGqS = ~V ■ (u • Vu) - C {t)afj \-u a \B + C { £\ , 
a a z \a p J 

S + 2H5- 4irGgS = -^V • (u • Vu) - \ [oV • (Su)}' + C {t)aP \-u a \ p + C* 



(17) 
(18) 
(19) 



Except for the presence of tensor- type perturbation, Eqs. (|17|) - (1T9|) are exactly the same as the ones known in the 
Newtonian theory. We note that these equations are valid in the presence of A. To the linear order, these are valid in 
the presence of general K, see Sec. lVIIll We have correctly identified the relativistic density and velocity perturbation 
variables which correspond to the Newtonian counterparts to the second order. In the relativistic context, our S and 
u are the density perturbation and (related to) the perturbed expansion scalar, respectively in the comoving gauge; 
the variables are cquivalently gauge-invariant. However, we were not able to identify relativistic variable which 
corresponds to the Newtonian gravitational potential to the second order; this is understandable if we consider the 
factor of two difference between Einstein's (post-Newtonian) and Newton's gravity theories in predicting the light 
bending under the gravitational field. To the linear order th e sp atial curvature perturbation in the zero-shear gauge 
can be identified as the perturbed Newtonian potential [HI, Il6j . Equations (fl"7| - (flT?)) include effects of gravitational 
waves to the density and velocity perturbations. Equations of the gravitational waves can be found in [l2j . 
To the third order, we have [l3| 



5 + -V • u = — V • (Su) + - \2tpu - V (A _1 X)1 • VS, 
a a a 

-V • ( u + -u ) + AnGfiS = — 3_V • (u • Vu) 



-^u.V(V.u) + lv. 



tp ( u • Vu — g u V • u 



-^[u.V(A-A) 



-u • VA 



3o ; 



(20) 



;AV-u{21) 



8 + 2^-5- AirGnS = [oV • (<5u)] + ^V • (u • Vu) + i {a [2<pu - V (A _1 A)] ■ V6} 



a 2 dt 

+ ^ u . V (V.u)-±V. 



tp ( u • Vu — gUV • u 



^ [u • V (A" 1 *)] - l u • VA - ^AV • u{22) 



where 



A = 2<pV ■ u - u • V^ + -A _1 V • [u • V (Vy>) + uAtp] . 



(23) 



In these equations we ignored the role of tensor-type perturbation; for a complete set of equations, see [13j. The 
variable ip is a perturbed-order metric (spatial curvature) variable in the comoving gauge condition, see later. 

All the third-order correction terms in Eqs. (|20 |) -(|22 p are simply of i^-order higher than the second-order rela- 
tivistic/Newtonian terms. Thus, the pure general relativistic effects are at least (/3-order higher than the relativis- 
tic/Newtonian ones in the second order equations. Thus, we only need the behavior of <p to the linear order which is 
related to the other hydrodynamic variables as 



-5$ + aA _1 V • u. 



(24) 



It also satisfies ITi 



= 0, (25) 
thus ip = C(x) with no decaying mode; this is true considering the presence of the cosmological constant, see (l7j . 
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IV. RELATIVISTIC FULLY NONLINEAR EQUATIONS 

In this section, for convenience, we present some additional covariant or ADM equations not available in [ill ]. In 
the multi-component situation, we have 



Tab = ^ ^(i)afc- (26) 



j 

The energy-momentum conservation gives 



T (i)a-b - hi)a, hi) a ~ °- ( 27 ) 



Tildes indicate the covariant quantities. 



A. Covariant equations 

We introduce the fluid quantities as 

T(i)ab = fj>(i)U(i) a U(i)b + P(i) (()ab + U(i) a W(i)f>) + q~{i) a U(i)b + q~(i) b U(i)a + ^(i)ab, (28) 

where 

u1i)U { ^ a = "I, Uf t) q il)a = = U( l) 7T (l)Qb , 7T (j;) ° = 0. (29) 

The fluid quantities of each component are based on the fluid four-vector Uu) a as 

= T (i)ab u (i)ul^, p{i) = -T^yabhfy, q^ a = -T{i)cd^i)\i)i> ^{i)ab = ^(i)cdfyi)a%)6 _ P(i)^(i)ab) (30) 
where h^ ab = g a b + u^ a n^ b . Equation ([27]) gives 

+ (m« +P(i)) %) + 9(i);a + q{i) h (i)a + ^(i)^(i)ab = -W(i)^(i) 0! ( 31 ) 
(A(i) +%)) + %)a + 9(i)b + %)6;c) + ^(i)ab + ^(i)a6 + g %) fyi)ab^ = \i)J{i)bi (32) 

where the kinematic quantities arc also based on the fluid four-vector u^ a as 

h(i) a h(i)bU(i)c;d = %)fafyi)b]"(i)c:d + %)(afyi)b)"(i)e;d = ^(j>6 + 9(i)ab = %)a:6 + ^(i)a^{i)b^ 

%) = u "i);a' ^(1)06 = <?(i)ab ~ g%)fyi)abi «(i)a = %)a;b**(i) = "(i)aj = A(i),a^fi)- (33) 

In the multi-component situation, we can derive the corresponding equation of Raychaudhury equation for the indi- 
vidual component. From u^a-bc — u^) a ;cb = U(i)dR d abc we can derive 

Ki) + - 5(i);a + " = 4^G (// - 3p - M^tifo^)) + A. (34) 

In the energy-frame we take q^) a = for each component of the fluids without losing any physical degree of freedom. 

In a single component situation, taking the energy-frame, the energy conservation equation, the momentum con- 
servation equation, and the Raychaudhury equation are [3] 

/»+ {fi + p)e + n ab d ab = 0, (35) 
(£ + P)5 a + /£(p, 6 + 7fg. c )=0, (36) 

+ ^ - ~a a , a + d ab a ab - Cj ab Cjab + InG (p + 3p)-A = 0. (37) 
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By combining Eqs. (j35|) - (|37|) we can derive 



A + K ah Qab 

A+A 



i / A + K ah d ab 

3 I fl + p 



4ttG (A + 3p) - A + a Qb a Qb - w afc £ Qb + 



ji + p 



(38) 



This equation was derived in Eq. (88) of [19( . see also [2(j. In the multi-component case, in the energy- frame, 
combining Eqs. (f3"Tj) - (f34f we can derive 



Hi) + + "(i) 1 ^)" \ 1 ( A(i) + *(^(i)ab + U a {l) I {l ) a 

3 



+ P(i) 



-4ttG (a - 3p - 2T afc U( l i) {t( j) ^ - A + a^a^ ab - u}fyCd( i}ab + 



A(») + P(<) 



(39) 



In [2l| Langlois and Vernizzi derived a simple covariant relation which leads to one of the conserved variable in the 
large-scale limit. These authors introduced 



C = h° 



Oih 



fl.b 



3(A + A). 

Using only the energy conservation in Eq. (|35p we can derive 



1 z 



3(A + A) 



P,a 



((M + p)9 



^A<» 





Ua^dbc 


7 (^),a 




[s(A + p). 


3(A + p) 



(40) 



(41) 



where £■& is a Lie derivative along u a with £u( a = (a-.bU b + ( b Ub-,a- Thus, for vanishing anisotropic pressure, we have 
the Langlois- Vernizzi relation [21| 



3(A + A) 



-^ha P,b ~ -A,6 



(42) 



/' 



These equations are valid in a single component fluid, or in multiple component fluids for the collective fluid variables. 
We can easily extend the relation to the individual fluid component as follows. We introduce 



C(i)a = tyi) 

Using only Eq. (|31j) we can derive 



a (i),b 



M(i),i> 



3(A(i) +P(i)) 



, a(j) = &(i^ a u°Q = 



(43) 



3(Aw +A(*)) 



P(i),a 



(Aw +A(i))^(, 



W(i)o (*(i) 5 '(i)6c+«(i)%)6 

3(A(i) +A(i)) 



3(Aw +A(i)) 



3(A« +A(i)) 



Thus, for vanishing anisotropic pressure and direct interactions among fluids, i.e., 7f(f) a & = = /(j) a , we have 



(44) 



P ? ti {i) I b ( - P(i) ~ \ 



(45) 



Application of these compact relations to large-scale conservation properties to the second order will be studied in 
Sec. [ 
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B. ADM equations 

The ADM formulation [22| is presented in Eqs. (2)-(13),(47),(48) of [HI- Interpretation of the ADM fluid quantities 
in Eqs. (45), (46) of [Tlj was based on the normal-frame fluid quantities; for relations to the energy-frame fluid 
quantities, see Eq. (|57[) below. The ADM fluid quantities of individual component are introduced as 



E W = h a n b f$ = A 2 T°°, J (i)Q ee -n b f (l) b a = NT (l) l 

S(i)aP = T(i) a p, S(i) EE h af3 S(i) a p, S(i) a p = S(i) a p — -h a pS(i). 

Equation ((2Tj) gives Eqs. (12), (13), (47), (48) in [ll|. Equations (10), (12), (47) in [11] can be arranged as 



(46) 



N 



(d Q - N a d a ) K% - - (K°) 2 = AttG (E + S) - A + K^K a(l 



N 



-N ]C 



K=(e+^ l(d -N a d a )E + ^(N 2 J a ). a -S^K a p 



K°t = E, 



(i) 



" 1 (9 - A^) E {1) + ± (a 2 J ( <* } ) a - S$K a p + i (j (i)0 - J (i)a JV«) 



(47) 
(48) 
(49) 



Momentum conservation equations for the collective and individual components can be found in Eqs. (13), (48) of 
[ill ]. By combining Eqs. ((17)) - ([4*9"]) we can derive the ADM counterpart of the density perturbation equations in Eqs. 



V. SECOND-ORDER PERTURBATIONS 

We use a metric convention in Eq. (49) of [ll[ 

goo ee -a 2 (1 + 2A) , g Qa = -a 2 B a , g a p = a 2 (g<® + 2C Q/3 ) . (50) 

The subindex indicates the conformal time 77 with adq = cdt. To the second-order in perturbation we introduce the 
fluid four- vector of individual component as 

fi W = ^0' 4) = - a ( 1 - A+ l A2 + \ V $) V V°> - B a V (t)a ) ; 



«(t)a = a (v (l)a -B a + AB a + 2V ( JC Q/3 ) , u (i)0 = —a\l + A — \^A 2 + ^V$V [i)a 



(51) 



In this definition of fluid four- vector we follow the notation in Eq. (53) of [Tl| . If we introduce u% = V^u^y we have 
Vfy = (1 — A)Vuy The fluid quantities of individual component are introduced as 

= + <fyi(*)j P(t) = P(i) + Sp(i), 7T( l ) Q /3 = a 2 Ii(i) a p, 7f(,) a = — a 2 II(j) Q/ g , 7T(i)00 = 0, (52) 

where from tt, .■. a = we have 

n w s - 2C^n (i)a/3 = 0. (53) 

The energy-momentum tensor of individual component in the energy-frame follows from Eq. (|28p as 

%)o = - <*M(i) - + *>(»)) - Bj V( l)Q , 

%° = +P(i)) ( V (i)<* ~ B »- AV {€)a + 2AB a + 2V^C a p) + (S m + Sp (i) ) {V (€)a - B a ) 

+n (i)Q/3 (vfa-B?), 

%s = (P(f) + %)) *f + (/*« + P«)) V% - B P ) + n (j) « - 2C^n (i)/37 . (54) 



9 



Using T b a = Tf~£, and the total fluid quantities in Eq. (82) of [Tl| we have 

A t = E^)' P = E p O')' 

j 3 

for the background order fluid quantities, and 

^ = E N) + faa +*&•)) - B ") - , 

+p) u Q = ^ +p (3) ) v (j > + (^(j) + ^o)) (Vy )Q - F a ) + n y £ (v (j)/3 - v^) 
n ^ = E { n w + (^(i) - Ba ) ( v w - v fi ) - \s$ 



(55) 



B^) v, 



U 



(56) 



for perturbed order fluid quantities to the second-order. From Eq. (|46|) we have the ADM fluid quantities based on 
the energy-frame fluid quantities 

E(i) = + S m + (V(?) - B ") ~ B ") > 

(MW +P(i)) - B a + AB a + 2V ( f ) C a/3 ) + + *P«)) (V(i)a - B a ) + n {i)a0 (v^ - B^j 



J(i)a = a 



S(i) = 3 (p w + 5p {i) ) + (HQ) +P(i)) (y$ - B a ) (V(i) Q - B a ) , 



<%)a/3 = a 2 S K(i)aP + (M(i) 



(V( 0a - B«) (V W/3 - Bp) - i 5 g - F") (V W7 - B 7 ) 



(57) 



We can compare Eq. l[57p with the ADM fluid quantities based on the normal-frame fluid quantities in Eq. (76) of 
[ll[; in a single component case we simply delete (i) subindices. 

In (llj the fluid quantities are based on the normal-frame vector. By taking Uu\ a = 0, uu\ a becomes the normal- 
frame vector h a , see Eq. (54) in fill ]. Based on the normal-frame vector, to the second order, the fluid quantities have 
contributions due to the frame choice: for example, even in the zero-pressure fluid, the perturbed pressure based on 
the normal-frame does not necessarily vanish to the second-order, see Eq. ([55)) below. By comparing Eq. ([57]) with 
Eq. (76) of [HI we have 

Srfi) = ^(f) + (MW + P«) (V(?j - £T) (V« )a - S a ) , 
SP& = + ^ (M W (vg - B Q ) (V (i)Q - B Q ) , 

Q(i)a = (At(0 - B« + AB a + 2V(JC a/s ) + + 5p w ) (V (i)a - B a ) + n (i)a/3 (v^ - B 13 ) , 

n") a0 = n (l)Q/3 + ( M(i) + p w ) (F (i)Q - B Q ) (F (i)/3 - B^) - \g% (Vfa - BP) (% 7 - B 7 ) . (58) 

For the total fluid quantities the relations between the two frames are presented in Eq. (87) of [1JJ]. Thus, by replacing 
all fluid quantities in Eqs. (99)-(107) of [ll[ using Eq. ([55]) and Eq. (87) of [ll| we have the equations in the energy 
frame. Using Qu\ a in Eq. (f58|) . Eq. ([54l) gives 



r (0 ° = (i-A)g Wl 



(59) 



As the fluid four- velocity of i-th component we can use either Qu\ a or V(j) Q — B a related by Eq. (p)£ 

The kinematic quantities in the energy-frame are presented in Eqs. (63)- (66) of [111 ]. In Eq. (|33[) we introduced 
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kinematic quantities for the individual component. To the second order we can show 



9 « = 4 - 4 A + \ c « + 5^1. + H A2 - 4 Bay ^ - \ Ac « + l v & + c ^ 



2 a 
2 



V(i)H/3) + Ca/3 - + (V(i)(a ~ #(a) (V (i )^) - Bfi) + V {l) ( a Aj i} + VfaC^ + 2V ( ] ) | (a C (/3)7 

a c«* (cr+^ )l7 ); 



0(i) a = Aa + - a (V(i)a - B a + AB a + 2V^C a/3 ) ' - 2AA, a - A- [a (V {i)a - B a )]' 







(0 



{V(i) a - B a ), + Bp\ a 



&(i)a0 = — V(f)0"(»)o!/3) CT(i)Q0 = 0; ^(i)aO = ~ ^)^(i)a/3, &(i)Q0 = 0; 5(»)0 = _ V(™)0(i) a , 



(60) 



where a prime indicates the time derivative based on 77. 

The gauge transformation properties of the fluid quantities are presented in Eqs. (232)-(235) of [ll[ for the normal- 
frame, and Eqs. (238) of [ll| for the energy-frame. A prescription to get the gauge transformation properties for 
individual fluid quantities is also presented below Eq. (235) of Under the gauge transformation we have 

*a w = sn (() - (V (i) + v w ) e° - s m , a e + \^U°e + (?e + re%) , 
sm = s P(t) - (V w + v w ) f + IpU^+p'o) (c°e ' + , 

V(i)a - B a + AB a + 2V^C aP = V [i)a -B a + AB a + 2V^C a p + £° Q ~ (V (l)a -B a )'?-^ (V (t)a - B a ) £° 

- (v m - Bfi ) e, a - (v {l)a - B a ) e + {a- e ~*°) e, a - ^% - ee\ a - 



n(i) a /3 — H(i)a/3 



n; 



(i)a/3 



2— n(i) a/ 3 ) £ - n (i ) Q/3:7 £ 7 - 2n (i ) 7 ( Q £ 1 



,0Y 



(61) 



A. Basic equations in the energy-frame 



The basic set of equations with fluid quantities based on the normal-frame is presented in Eqs. (99)-(107) of [Tll | . 
By using Eq. (]58[) we can recover the equations with fluid quantities based on the energy frame. For convenience, in 
the following we present the complete set of equations with fluid quantities in the energy frame. These equations are 
written without taking any gauge conditions yet, thus in a sort of gauge-ready form. To the linear order this method 
was suggested by Bardeen |23l |24|. 
Definition of 5K: 



K + 3H + 5K- 3HA + + -B a v 



= -A [ -HA - CI 



-B c 



| a 



HB a B a + ±-B a (2Cf |/3 - G| |a ) + 2C a/} (d aP + ^B aW 



no. (62) 
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Energy constraint equation: 



WirGfi + 2A- 6H 2 - + lGnGSfi + 4H5K - \ f 2Cf p - 2Cf% - \r {z) C l c 

Cl Qi \ o 



\5K 2 - 16ttG (h + p) (V a - B a ) (V a - B a ) - C a 



-B( a \f3) 



C afj + -B al 



Cl + -B c 
a 



+ — 



= m. 



(~Cl m - Cl hfJ + Cj\ + C^) + -R^C«Cl 

(2Cj h C^) (2C af \ a - C^) + C^ (3C aPh - 2C aiW ) 
Momentum constraint equation: 

= ~ASK, a - 8nGa [(Sf* + Sp) (V a - B a ) + (fi + p)A (V a - B a ) + U a p (V p - B (3 )} 



(63) 



1/3 



B- 



|/3 3 ^ V 7 a It 



Trace of the ADM propagation equation: 



la 



|a 



= "-2c 



(64) 



A 



3ff + 3iJ 2 + AttG (p + 3p) - Aj + SK + 2^(5^ - 47rG (fy + 36p) + [ 3H + — ) A 

= A5K - -5K, a B a + Uk 2 + 8ttG (jji + p) (V a - B a ) (V a - B a ) 
ci o 



\H (3A* - B«B a ) + 1 



2AAA + A- a A, a - 1a i K B a B a ) + A a (2Cf |/3 - C| |a ) + 2C^A , 



- - + --B' 



(65) 
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Tracefree ADM propagation equation: 
1 



3H 



-A' 



— r\ a 

3^ 



1 



+3H 



+3 ^(^ + H)-> 17 



87rGn° 



A + 2C ( C, 7 + lB Wl) ) + - a B, (Cp + C^ a - <fp 



A + 2C'"[C i - + - B{ph) )+±B y ( 



° |/3 "+" °/3 °/3 



B~< + SK 



+ - 



+3iJ 



^4'% + \ {-A 2 + W% - 2CAj, h (Cp + 0$ a - Cf<) A, 



C 7 + ^B 7 , 7 ] A + 2C^ S (c 7S + -B, 



- ( C 7 + V, 7 ) A - i (&< + ^ + SK + i n 



1 



-AAl^ 7 + l - (-A 2 + B S B S ) h 7 - 2CT^, 7 | 5 - (2C^ |7 - Cf) A, 5 



a^l V 51/37 5,3 7 ~~ /3|<57 ~ 57 0j +ZU V°7l/3<5 + U 0hS ~ °/3 7 5 ~ °<5| 7/ 3 

+ (2C 7 |7 - C 7 5 ) (C as |/3 + - C^) - C^C 7 ^ + 2C« 7 I 5 (C 0Sh ~ C 01{5 ) 

-±s% 4C^(c* lSe + c; ]€5 -c^\-ct h5 ) -±rWc$C2 

+ ( 2 Q |e - C| |5 ) (2C 75 |7 - Cf) + (2C\ els - 3C 7 



3 



-16nGC a m Pl + 8irG (p + p) 
Energy conservation equation: 



(V a - B a ) (Vp B P ) - -Jp (V B 1 ) (V, B 7 ) 



Hp- 



(66) 



[A + 3£T (/x + p)] + <5A + 3i/ (J/x + dp) -(p+p) (SK - 3HA) + - + p) [V a - B a + AB a + 2V 13 C^} |q 
= - U[i. a B a + {5n + 5p) (SK - 3HA) + (ji+ P ) ASK + '^H(n+ p) (A 2 - B a B a ) 
~ [a 4 (p+p) (r - B a ) (V a - B a )} ■ - 1 [{6n + Sp) (V a - B a ) + (V fj - Bp)} |q 
+ \ ( M + p) {-A (V a B a ) [a 2A, a (V a -B a ) + 2 [C^ (Vp - Bp)] ]a - C^p - B?) } 
-11"? (d a p + \B a ^ = n 5 . (67) 



13 



Momentum conservation equation: 

\ [a* (M + P) (V a -B a + AB a + 2V> 3 C a p)] ' + - (m +p) A. a + - (s P , a + Ilf , 

Gj Gj Ob ^ 

= (M + P) (SK - 3HA) (V a ~B a )-^ {a 4 [(5[i + 5p) (V a - B a ) + ng (Vp - Bp)]}' 



+- 1 - (Sp, a + IF |/3 J A- {6(1 + Sp) A a 

- (p + p) [-AA, a + Bp la V? + (V a - B a )\p V 13 + (V a - B a ) (V p - B 13 )^ 

+2 (c^n Q7 ) |/3 - c^w a + c^X - a il^ = n 6a . 



(68) 



In the multi-component situation we additionally have the energy and the momentum conservation of individual 
component. Using the energy-frame fluid quantities, Eqs. (106), (107) in [ll| become 

£(») + 3# + p(i)) + - J(i) + ^A(i) + 3B (^(i) + *P(i)) - (M(i) + P(i)) (Mf - 3BA) 



+- (mw + p«) 



Vg -B« + AB a + 2U ( f ) C^ 
I . . , , .... 3 



la a 



(i)0 



^ (l) , Q B Q + (<5 M(l) + 5 P(l) ) (5K - 3BA) + +p (<) ) A&K" + -B ( M(<) + (A 2 - B a B a ) 
a 4 (/i(i) + P (i) ) (v% - B Q ) (V (i)a - B a )] ' - \ + (^W - + < - 

+ ^ (AW + P«) {-A (V(?) - B°) |q - 2A, Q (vft - B«) + 2 [C°"» (V (i)/3 - B,)] |a - (V (i)/3 - B^) 



4 + (V (i)a - B Q + AB a + 2Vf i) C a0 



(69) 



+ - +P(i)) Aa + - (%),a + n ( 4 )a| /3 - M (i)a 

= (m« + P«) (MT - 3BA) (V Wa - Ba) - ^ {a 4 + tfp (<) ) (F (i)a - B a ) + n (i) f (V (i)/3 - B^)] }' 



,a + Bflat$ + (V Wa " B«) ]ff + (V (l)a ~ B a ) (v* - B?) 



+2 (c^n (l)ai ) - c| |7 n (i) 2 + c^| a n w ; - A/?n (i) f = n (l)6c 



(70) 



By removing indices indicating the components in Eqs. (|69j) .(|70 |) we recover equations for the collective component 
which coincide with the equations in a single component situation in Eqs. (|67p.(f" 
To the background order, from Eqs. (|6"9"|) , (|6"5)) we have 



(t(i) + 3H (mq +P(i)) = — I(i)o, 
ft + 3H(n+p) = 0, 

3B + 3B 2 = + 3 PU)) + Ac2 > 



8nG 



Kc 2 Ac 2 



m - — + — 

3 



(71) 
(72) 
(73) 

(74) 



where we recovered the speed of light c. Dimensions are 

[Gq] = T- 2 , [c] - LT~\ [rj\ = 1, [p] = \p] = [gc 2 ], [a] = L, [K] = 1, [A] = L~ 2 , [7 (i)0 ] - [//]. (75) 
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Equation (|72f follows from the sum of Eq. (|7T|) over components. Equation ([74]) follows from integrating Eq. (|73|) 
where i-T-term can be regarded as an integration constant; in Einstein's gravity K-term can be normalized as the 
sign of spatial curvature. Compared with the Newtonian background equations in Eq. ([7|), ignoring the direct inter- 
action terms in Eq. (|7Tj) . the presence of pressure terms in Eqs. (|7ip - (jT4")l is the pure general relativistic effect. The 
cosmological constant A can be introduced by hand even in the Newtonian case. 



B. Decomposition 

We decompose the metric to three perturbation types 



A = a, B a = f3. a + flW , C aP = <pg®> + 7,^ + Cgj^ + C7<g, 



(76) 



where superscripts (v) and (t) indicate the transverse vector-type, and transverse-tracefree tensor-type perturbations, 
respectively. We introduce 



X = a (0 + c^aj) , 9™ = + c^aC™, 
which are spatially gauge-invariant to the linear order. We set 

K" = -3H + k. 



(77) 



(78) 



We will identify n with Newtonian velocity variable which will be an important step in our analysis, see Eqs. 
For the fluid quantities we decompose 



\i) a = a f V(,-) Q — B a + AB a + 2Vf^C a p) = av {i)a = a + 



n 



(i)a0 = -5 



H-(i),a\/3 ~ o9al AU (i) ) + - U \i)( a \/3) + n (»ja/3' SI (i)<* = 5I (i),<* + SI (i)a 



A") 



(*) 



(79) 
(80) 



The perturbed fluid velocity variables and subtly differ from the ones introduced in [TTj]; see Sec. IV Cl For 
the collective fluid component or for a single component case, we simply delete (z) subindices. For isotropic pressure 
we introduce 



5p = c 2 s 5[i + e, c 2 s = -r- 



(81) 



The perturbation variable e is called an entropic perturbation. Defined in this way e is gauge-invariant only to the 
linear order. To the second order, from Eq. (|6Tj) we can derive the following gauge-invariant combination 



_ 6(i 

opsu = e - 

A* 



(82) 



In our notation, (^,5^ is a gauge-invariant combination which is the same as dp in the <5/i = slicing (temporal gauge) 
condition to the second order; as the spatial gauge we take 7 = to the second order; for the derivation, see the 
prescription below Eq. (266) of pjjj. In the multi-component case, we similarly have 



, _ <*M(*) 

°P(fiS m = e(i) - -r 



(0 + ^(4)) S V(i 



where 



We have \2l 



a — 2 x 1 2 _ 

S P(i) = c (i) ^ ( i) + e (i) , c (i) = — . 



e = e rcl + e int , e re i = (cf^ — cf) <fyt(j), e int = ^ e 
i 



0')- 



(83) 



(84) 



(85) 
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Equation (f56|) gives 

si* = [ Sfl ti) + (m +p(j)) v u) ( v ti)<* - u «) 



<^oo + 3 0*y) +PU)) v u) («o> - v ») 



3 



(86) 



By recovering c, dimensions of the variables are 



[g«»] = [g ab ] = [u a ] = i, [T ab ] = H = i, [v] = [A] = i, M = [c 2 ] = l, 

[A] = [B Q ] = [C Q/3 ] = 1, [a] = [p] = [/?] = [ 7 ] = [BW] = [CW] = = [<7<g] = 1, [ X ] = L, 

= N = [e] = [n a/3 ] = [u%] = [<5] = [v^] = i, M = [vW] = i, [n] = L 2 [ M ], [nW] 

Scalar- type perturbation equations can be derived from Eqs. (f6"2"| - (170[) 

ft — 3Ha + 3ip H — »v = no, 
or 

4 „ A + 3if 1 

47rGo/.i + ilK H ^ — w = —rii, 

a z 4 

A -|- 3_Zi" 3 

K H 5 — x - 127rG f/i + p) av = n 2 = - A -1 V a n 2a , 

er 2 

K + 2Hk - AttG (Sfi + 3Sp) + \3H + a = n 3 , 

X + Hx-^-a- 8nGU = n 4 = ^a 2 (A + 3Ky l A -1 V V /J n 4 f , 

A 



A- 1 V a n 6c 



■-L\x\. (87) 



+ 3H (Sfj, + 8p) — (fi + p) k — 3Ha H v I = n^, 



[a 4 (fj, + p)v\ 1 1 



a 4 (fx + p) a a (/x + p) 



2A + 3if TT . 
*p+-_^II]=n 6 



M + p 



+ 3 # (^M(i) + <5p(i)) ^ +P(i)) ~ % Ha + ~~ V (*^J + ^%)° = n5 M 
[a 4 (p (l) +p w ) « W ]' 1 J 2 

— V \ d P(i)+o 



2 A + 3AT 



= s ~ - * A~ 1 V Q! n 6o ,(i). 
Equations for the vector-type perturbation follow from Eqs. (|6^ . (l6l^ . (|68p . (fT0")) 

A 1K + 8irG(n + p) V M = - (n 2a - V a A^ V^ 2/3 ) ee , 



2a 2 

^ } + 2ff«M - 87TGIIW = 2a (A + 2/f)" 1 (v^n 4 f - V^^V^V 
[o*(m+p)w£ o) ]- A + 2^n^ } 



(89) 
(90) 
(91) 
(92) 
(93) 

(94) 
(95) 

(96) 
(97) 



„ 0\- „(") 



1 



a 4 (/x+p) 2a 2 H+p H+P 

-Af,, i ^. , \„S V ) 1- « . „ T ^ tt( 1 ') 



- («) 



n ea - V Q A 1 V l3 n 6f3 ) = ,> 6a , 
[o 4 (/*W+P(i))«(i)a]' A + 2K n^ a i (JjgJ, 



a 4 (^) +i>(i)) 



2a 2 V(i)+P{i) 



(n 6 (») a - V a A ^n^, 



)/») = n 6( 



(99) 



Uioo) 
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Equations for the tensor-type perturbation follow from Eq. (|66[) 

r-w + wrW A-2K w (t) 



-2V (Q (A + 2#) 1 (V^n 4/3)7 - V^)A- x V^Van4) = n^. (101) 
In order to derive eqs. (|92I98I101|) it is convenient to show 

^ (Va V, - A ) + SttGII) + 1 (a 2 *[i)) ' - 8^x1^ 

+Cg + - ^^C% - 8nGH% = n 4a p, (102) 

which follows from Eq. (f6"6"| . Quadratic combinations of linear-order perturbation variables of all three types of 
perturbations contribute to all three types of perturbation to the second order. 



C. Comoving gauge and irrotational condition 

In Eq. (180) of jTTj] we introduced 

Q(i)a = +P(i)) (-%,<* + i ( 103 ) 

where we put overbars to and v^ a in order to distinguish these from our new notations to be used in this paper; 
in [ll[ we didn't have overbars. From Eq. (|55|) we have 

-v {i) , a + B w ^ - s Q + ab q + 2^c a/3 + + fr« {v _ Ba) Jk^ fy& _ B p\ (104) 

It is more convenient to introduce the decomposition in Eq. (|79[) . Thus, we have 

«W,a + «(i)a " »W.« + + — ^ «W,a + « (i) aJ + — T~ ^ «(() + J • (^J 

The variable Vu) a introduced in fllj cannot be regarded as a proper vector-type perturbation. However, if we also 
take the temporal comoving gauge in [ll[ which sets v = together with u£\ Q = 0j we have v = = u£\ a ! these 
are the same as taking the proper irrotational condition (f/X a = 0) and the temporal comoving gauge (« = 0). Our 
analyses in [H|, [l3| are, in fact, based on taking these two conditions together. In the irrotational fluids, the temporal 
comoving gauge v = leads to u a = 0, thus u a coincides with the normal frame four- vector h a - 
From Eq. (|6H|) we have 



(v) a ( _ (v) 



2A + 3X \ A + 3K (u) 

(106) 



where we used Eq. (|70[) to the linear order. For vanishing vector-type perturbation we set v$ a = 0, etc. In this case, 



we have 



,,-.„< = -7r-7T-%[a (%) + l^-T^ n W -%)) » ( 107 ) 

■ V a a j 3] 



V a 7,/3] 

which vanishes for 8pu\ = = ILj) and = 0. 
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VI. EFFECTS OF PRESSURE 



We consider a single component situation without rotational perturbations. Equations (|65p . (f67|) . ([68|) , and Eq. (|64l) 
to the linear order provide a complete set of equations we need in the following. 



A. Irrotational case 

As an irrotational fluid we ignore all vector-type perturbations, thus i4 — — Ca = Ha = 0. Quadratic 
combinations of linear-order perturbations of all three-types of perturbations contribute to each type of second-order 
perturbations. Thus, concerning the second-order scalar-type and vector-type perturbations, ignoring the pure vector- 
type perturbation corresponds to ignoring the vector-type perturbation only to the linear order; to the linear order the 
vector-type perturbation only has a decaying solution in the expanding phase. Assuming the background equations 
are valid Eqs. §7$ , |JB5J| , (JBSJ) , JM|) give 



S + 3H ( ^ - ,rS 



(1 + w) ( k H v - 3Ha 



~S,a(3' a + (s + - 3Ha) + (1 + to) \k + ^-v ] 



+\H (1 + w) (a 2 - (3>°<f3, a ) - (± Xia]fl + ~l(l + w) ( w g + 7 , q|/3 + C«) 



+- (1 + iu) (a [H+-^- 
a IV M+P 

+ — \(Su + dp) w a + Il a Pv g] 

a/i L ' ' J I' 



(108) 



- i [a i (fx+p)v ja ] -- (fi+p)a, a + Sp, a +n^, g = -((i+p)(v - P)' 13 (v- 0) 
a* a L "i p J a ' 

+ (fi + p) v j0l (k + —v - 3Ha) - i {a 4 [(tfyi + 5p) u Q + EE^u^] } 



1 



/3 



(fyx + <5p) a a + a (n+p) a, Q — 5p, a - ^ a 



1/3 



+ 2 (c^n a7 ) - n^, + n^ca J , (109) 



« + 2Hk - AnG (5(i + 35p) + ( 3H + 4 I a = --« c# a + «/i + ^k 2 + ^i? (3a 2 - /3' a /3 Q ) + 8ttG (u + p)v a v a 

\ a z J a 3 2 

+ ^ |2Aa (a + ^) + a> a [a - ip + (A + 4ftT) 7 ]_ Q + 2a'«l' 3 (7^ + C<g) - ^A (/?' Q /3 Q ) J 

>* + e< ' W )(^.* + cS)-i(^) 2 , (110) 



A + 3if 



X + 127rGa (/x+p)v. 



(Ill) 



Equation (jllip is valid to the linear order. 



B. Comoving gauge 

As the temporal and spatial gauge conditions we set 

v = = 7, 

thus, j3 — x/ a - Under these gauge conditions Eq. (|109p gives 

1 „ c? 



2a z 1 + w 



2(1 + tu) . 



an, 



(112) 



(113) 
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where the imperfect fluid contributions (stresses) are denoted by an with 



1 



P+P 
1 

p + p 



2 A + 3A' 

e + = — n 



3 a 2 



, 2 c 2 X 
-1 + - — —8- 



1+W pk+P 



2 A + 3X 
e H - — II 



A _1 V" 



l + w 



3 a 2 
H+p 



iii + 2 (c^n Q7 ) |/3 - n^cf |7 + n^cg |a 

(114) 



Using this, Eqs. (fT08| - (fTTTT) give 

1 3 r 2 

S - 3HwS -(1+w)k= — nS,aX' a + ^ + ~H-^5 2 + S n 



k + 2Hk - 4ttG^S 



cj A + 3K 
1 + w a 2 



2 1 + k; 



5 = 5" K ,aX' Q + o K " 



(115) 



2(1 + 



4ttG^ 



1 + 2c 2 A + 3K 
1 + w a 2 



1 + w 



3 (a 2: 



-«n, 
A + 3X 



(116) 
(117) 



where Eq. (|117p is valid only to the linear order, and 



<5 n = 2H 



A + 3K II e 



-[ k+3H 



1 



Ai + P 
1 



a' fj, fj, 

2 A + 3X 
e H = — II 



l + w 



3 a 2 



3 „1 



3H5 - (1 + w) k + -H- e 
2 ji 



2 A + 3K 

3 ? 



n 



-3fT-A _x V a 
Kn = \2nGe ( 1 - 



p. + p 



nf + 2 (c^n a7 ) - n^ |7 + u^c; [a 



l + w 



! 3ff +4 



1 



a J M + P 



2 A + 31T 
e+ 3^^ n 



-1 



2c 2 



l + w + P 



2 A + 3X 
e H — II 



3 a 2 



c 2 , ^ + n^ |7 



m + p 

2(c^n a7 ) |/3 -n^c| |7 + n^ 

2Fk + AttG (1 + 6c 2 ) ^ + 127rGe 



-5 



l + w 



>/3 



P + P 



t7 



c 2 3K „ 3- 1 
/z - 



2 A + 3if 



1 



1 1 

/ilpo 2 



2 A + 3X 
e H = — II 



3 a 2 



1 A 

fi + p a 2 



l + w a 2 2 H + p 

2A + 3/0 ri I 

e H - — 11 



3 a 2 



2^, + p 



3 a 2 

2A + 3K 



3 a 2 



n ) +2ip- 



n 



A 



l + w 



2 A + 3A, 

3 ^ 



n i/j a - 2 e + - 



2 A + 3K 



3 a 2 



n 



(118) 
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Combining Eqs. (|115p . (|116[) we can derive 
1 + w 



a 2 H 



H 2 fa 3 fi 



+ - 



(/i + p)a V H 
1 + w 



A 



1 ( 2„ x X ,.,a i 3 „2tr c s ^2 



1 + w \ 2 1+w 



1 



1 + 2c 2 A + 3K 



+c 2 ^ (2a 2 ^ K + ^ >a - 2^A5 - 2S' a WC<$) +(l + w) C (t)a/3 ( JfX,a|/3 + C% 



ft ^ 1 + w / a 2 r , 
+ (1 + to) Ku H 5— -r- — % o<>n,aX' 



(119) 



Notice that the equations above are valid in the presence of general K and A. 



C. Newtonian correspondence 



Guided by our success in the zero-pressure case, we continue to identify 



-V • u, 

a 



(120) 



to the second order. Using Eq. (jTTTJ) , assuming K = 0, we can identify 



A 1„ c ^ 

K = -c-sX= Vu ; u=-Vx = -cVv x , 

er a a 

to the linear order. We have recovered the speed of light c. We have [u] = LT^ 1 . Equation (|113[1 becomes 



(121) 



1 r 2 

2c 2 1 + w 



1-4+^,5 
2(1 + w) 



an- 



(122) 



Equations (flT5|) . (fTT6|) give 



1 1 3 r 2 

5 - 3wHS + (1 + to) -V • u = — V • (<Ju) + --^H5 2 + S n , 
a a 21 + 10 

-V-(u + ffn) 

a 



S + ^-c 2 4<5 = 4v ' (u • Vu) - ( - 

1 + w a a 2 \a 



1 + w 12 



1 / AnG^i l + 2c 2 2 A 



+ 



1 + to a 



c 2 ^ S z + 2HS-V ■ u + -=■ 



2 V AS - (V(f) ■ VS + 2S' a ^C ( *l] \ - k4L24) 



(123) 



We have [an] = 1, [Sn] —T 1 , and [«n] =T 2 . Combining these equations or Eq. (|119[1 gives 



1 + w 
a 2 H 



H 2 



a 3 /j 



(M + P) a V H 
1 + w 



2 2 A A 

C S C -rd 

a 2 



a 



, V(«.Vu)-i+^ 



1 + to 



2 

s c2 



v ; 2 1+w 



}' + (l + -)^(^ Q „ + ^) 



1 2 /4^Gm _ 1 + 2C 2 . 2 A\ _ 

2 8 ^ c 2 1 + w aV s a 2 



2ai2W • u + 2^c 2 A<S - c 2 (Vy>) • + 2c 2 S' a ^C { ^ 



+ (1 + w)kh + 



1 + w / a 



1 



Sn] u V5n- 

1 + w J a 



(125) 



We note that, to the linear order, Eq. (|125|) is valid in the presence of general K and A, see Eq. (|119[) . 
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D. Linear-order relativistic pressure corrections 

To the linear order, ignoring the entropic perturbation e and the anisotropic stress II, Eqs. (|123[) - 1|125|) become 



S - 3wHS + (1 + w) -V • u = 0, 
a 

-V ■ (u + Hu) + AttGqS = 4 — > 

a t + w a z g 



1 + w 
a 2 H 



H 1 



a(l + w)g \ H 



a 3 g 



A Sp 

9 ' 



(126) 
(127) 

(128) 



where ignoring the specific internal energy density e we used /i = gc 2 , thus 5 = <5f?/f?; in general we have /i = g (c 2 + e) 
[l8l |. Equation (|128p can be expanded to give 



1 + w 
~aJW 



H 2 



a 3 g 



a(l + w)g \ H 



(2 - 6w + 3c 2 ) HS 



A Sp 

9 

a g 

(l + 8it;- 6c 2 - 3w 2 ) AnGg - 12 (w - c 2 ) 



Kc 2 



(5w-3c1)Ac 2 



5, (129) 



which is valid in the presence of K and A. Equation of density perturbation in the comoving gauge was first derived 
by Nariai in [H]. 

In a single component case the Newtonian equations in Eqs. (|8"l), (fT2"l) . (116[) give 

(130) 
(131) 

(132) 



5+ -V • u = 0, 

a 

1 A Sp 

-V • (u + Hu) + AnGgS = 7 — , 

a a 2 g 



5 + 2H5 - AnGgS = 



1 



a 2 H 



a 2 H 2 



H 



A Sp 
a 2 g 



Comparing Eqs. (|126[) . (|12T[) with Eqs. (|130p . (|13ip . we notice the presence of w = p/{gc 2 ) term in three places in 
the relativistic equations. Even to the linear order the presence of these terms should be regarded as pure general 
relativistic effect of the isotropic pressure. The effects of pressures to the second order compared with the Newtonian 
equations can be found in Eqs. (|123[) - (|125p which should be compared with Newtonian equations in Eqs. (|8")). (fT2"l) . (116p . 
Pressure has the genuine relativistic role in cosmology even in the background level. 



VII. EFFECTS OF MULTI-COMPONENT 



We assume zero-pressure medium, thus set 



P(i) = 0, = = U(i) a p. 



(133) 



From Eq. (|86p we notice that to the second order the collective fluid quantities differ from simple sum of individual 
fluid quantities. Even in the zero-pressure mediums, we have 

si* = [^(j) + m v u) (%> - . 

3 

j 
3 

U P = J2 ^CJ) V (3) ( V i3)t3 ~ V 0) ~ o S P V (j) ( V (3h - V ~f) . ( 134 ) 



thus Sp 7^ Ug to the second order. 
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Equations @ , JTOJ) , lESJ , (E3) , JSH1> , lEl give 

1 



>(i) 



SK + 3HA + - 



1 + «gJ = + (MT - 3HA) + ASK + ~H (A 2 - B a B a ) 



■«Wa (2*3) + ~ ^)|a - + ^ (C°%)/j)| a ~ 



i [a (1 + u Wce ] ' + -A, a = (SK - 3HA) v {i)c 



AA. a - S {i) A, a - B p B (j \ a - (w(i)aujflj - w (i)Q | /3 B /3 - v^B ^ 



5if + 2iJ<^ - 4ttG (5u + 3M + [3H + ^]a = ASK - -SK a B a + -SK 2 + -H (3A 2 - B a B a 

\ a z J a 3 2 

2AAA + A> a A >a - (B a B a ) + (2Cf |/3 - Cg a ) + 2C^A, a | / 
+ (c aP + -B^) (c a/3 + -B al s) - ~ (ci + -B a la ) + 87rG^v a v a 



a *W] 3 



S-SK + 3HA + - [(1 + 5) v a ] la = —S, a B a + 5 (SK - 3HA) + ASK + -H (A 2 - B a B a ) 

-v a (2v a + Hv a ) - -Av a {a - -A, a v a + - (C af3 v p ) - -C^% - 3H^, 

- [a (1 + S) v a V + -A a = (SK - 3HA) v a 
a a ' 



AA <a - S l{) A, a - B fi B p \ a - (v a v ) - v a{0 B - vpB - - (s P , a + II 



+ B W 



l(c2 + -B 1 



If 



h 



2 x 

-8K, a = -8-KGa } ,tlU)V(j)g- 



where Eq. (|140|) is valid to the linear order. 



(135) 



(136) 



(137) 

(138) 

(139) 
(140) 



A. Irrotational case 



Assuming an irrotational condition we ignore all vector-type perturbations. As the spatial gauge condition we take 

7 = 0, (141) 
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thus, f3 = x/ a - Equations (|135|) - (|140[) become 

1 



Hi) 



K + 3Ha 



(1 + 5 {l) ) v$ 



> s (i), a r a + (« - 3Ha) + aK+-H 



11 2 2 



~ [a(l + <yw)«(<)]; a - ~a, a 
= (« - 3fla) u/i) Q + - 



1/3' 



k + 2Hn — AttG (Six + 3Sp) + [ 3H H — - ] a = ak -k a \ 



1 2,3 • / 2 1 



1 



2aAa + a ,Q a , 



A 

2^2 



(x ,Q X,a) " + 2(^Aa + 2C« Q/3 a, 



>\0 



+ fa'** + 3X' a|/3 ) fa-8 + ixaU,) - J 2 + Q 



3 \a 

S-k + 3Ha - - [(1 + S) v a ]. = -^5, a x a +S(k- 3Ha) + a K +^H (a 2 - \ X ' a X, a 
a a z 2 \ a z 

+Hv a v a + -aAv + -cp' a v a - -ipAv - -C (t)af3 v a \ - 3ff — , 
a a : a a /i 

— [a (1 + 8) v] a a. a = (k — 3Ha) w a 



1 

+ - 

a 

A + 3K 



X + K = 127rGa^^ (j) w (j) . 



Equation (I147|) is valid to the linear order. 



(142) 



(143) 



(144) 

(145) 

(146) 
(147) 



B. Linear perturbations 



To the linear order Eqs. (|142|) - (|146[) give 



<5(ij — k + 3Ha — c — = 0, 
c 

V(i) + Hv (i) a = 0, 



A 



k + 2Hn - 4vrG^ Q(j) s U) + 3# + c 2 — a = 0, 



A + 3if 12ttG v-^ 

c s — x + k= __a^ecfl«W)> 



A 



<5 — K + 3Ha — c — v — 0, 
a 

c 

u + -ffw a = 0, 

a 



(148) 
(149) 

(150) 

(151) 

(152) 
(153) 
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where we have recovered the speed of light c. The following additional equations can be found in Eqs. (195), (196), (199) 

of mi 

k - 3Ha + 3<p + c4-X = 0, (154) 
a z 

A-kGqS + Hk + c 2 + ? K ip = 0, (155) 
or 

-(x + Hx)-<p-a = 0. (156) 
c 



Equation (|156[) gives 
Equation (|153[) gives 
Equations (fT5T|) , (fl"54|) give 

Thus, for K = we have 



■x = -fx- ( 157 ) 
a v = 0. (158) 



Kc 

<Pv = —oXv (159) 



= 0, (160) 
which is valid even in the presence of multi-components and the cosmological constant. Equations (|149f) . (|153[) give 

i> x + Hv x + ~<p x = 0, (161) 
c 

V (i)X + Hv (i)x + -<Px = °' ( 162 ) 

Equations (fT5l|) , (fT53)) give 

c 2 ^_3K ^ = _ AitGq5v = -AitG^q {j) 8 {j)v = -4ttG^ ■ (163) 

j 3 

We can derive density perturbation equation in many different temporal gauge (hypersurface) conditions all of which 
naturally correspond to gauge- invariant variables. In a single component case, density perturbation in the comoving 
gauge (v = 0) is known to give Newtonian result. In the multi-component situation we have many different comoving 
gauge conditions. Here, we consider two such gauges for 5^ variable: one based on v — gauge, and the other based 
on vm = gauge for a specific I. 

(I) Equation (|148[) evaluated in the v — gauge, and using Eq. (|151[) we can derive 

• A 3K . , 

0(i)v ~ c — v (i)x - C V X = °) ( 164 ) 

& a 

where we used \v = X ~ av = — a v x , v V(i) = v — v^, and a v = a — c _1 (av) — 0. Using Eqs. (|16ip . (|162p . (|163p . we 
have 

8(i) v + 2H5 {l)v - 4ttG^ Q{j)5{j)v = 0- (165) 

3 

This coincides exactly with the Newtonian result in Eq. (flUI) to the linear order, even in the presence of K. Thus, we 
may identify Su\ v as the Newtonian density perturbation Si to the linear order even in the presence of K. For K = 
we may also identify -cVi)(,) x as the Newtonian velocity perturbation u^. However, in the presence of K, we cannot 
identify the relativistic variables which correspond to the Newtonian velocity perturbation of individual component. 
Therefore, to the linear order we have the following Newtonian correspondences 



Si = u t = -cVu(i) x , 



(166) 
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where the latter one is valid only for K = 0. 

(II) Evaluating Eq. (|148[) in the — gauge for a specific £, and using Eq. (|I5I[) we can derive 

• A + 3K I2ttG , , A , s f10tT . 

*W»c« _ c a = —^- a 2^ey) {v(j) - v w ) + c— - , (167) 

3 

where we used Xv w = X ~ av (e) = - av {l)xi «(<)« w = - and = a - c" 1 (av^)) = 0. Using Eqs. 

([I6l l) .(p2 l) .(p3 l) . we have 

S (i)v(e) + 2H5 (l)v{e) - 4ttG^ Q(j)8(j)v W = — ^— a# (v (e) - v (f) ) . (168) 

3 3 

The terms in right-hand-sides of Eqs. (|167ft . (1168ft look like relativistic correction terms present even to the linear 
order based on the variable S^ v . Since no such correction terms appear in Eq. (|165[) based on the variable S^ v , 
the relativistic correction terms in Eq. (|168|) can be regarded as being caused by a complicated hypersurface (gauge) 
choice. 



1. Exact solutions 



Assuming K = 0, we can identify Newtonian perturbation variables as 
„ „ 1. 



-Vu, u 



-c f x , 



-cVw 



(Ox- 



Equations (fl^ . (fl^ .(fl ^ . (fT4^ . (fT63l) become 



5 = — Vu, 

a 

u + Hu= --V6&, 



Si = — V • u i7 
a 

\ii + Hui = — V<5$, 



\5<S> = AnGgS = iirG^QjSj 



(169) 

(170) 
(171) 
(172) 
(173) 
(174) 



Under the identification in Eq. (|169[) these equations are valid in both Newton's and Einstein's gravity theories. 
Equations (fT7D|) . (fTTTj) , fl7lft . and Eqs. {T72J| , {T73JI , (HZU) , respectively, give 



S + 2HS - AnGgS = 



1 



a 2 H 



a 2 H 2 [ - 
H 



= 0, 



Si + 2HS t - AttG Qj$j = 0. 

3 

Equation (| 175ft has an exact solution 



5(x,t) = H 



Cg(x) 



+ c d (x) 



(175) 
(176) 

(177) 



where c g and Cd are integration constants which indicate the relatively growing and decaying solutions in expanding 
phase; we do not consider the lower bound of integration which is absorbed to the Cd mode. Equations (|175[) . (|176[) . 
and the solution in Eq. (|177[) are valid considering general K and A in the background world model. Equation (1 174ft 
can be solved to give 



8>S> 



-Gga* 



(178) 
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From Eqs. (fTTOf . (fTTTj) . (fT74| we can show [3 

/ VJ$ Y 1 

u=-a( i ^J+-D W) V-D.O, (179) 

where the D term is the solution of the homogeneous part of Eq. (|17ip ; it decouples from the density inhomogeneity 
and corresponds to the peculiar velocity in the background world model. Since the D term is not connected to the 
density inhomogeneity and simply decays, we may ignore it to the linear order. 
Now, for the individual component, from Eqs. (|172jl . ()173|) we have 

u. t = -a ( 2 ) +-Vd 4 (x) + -D i (x), V-Di = 0, (180) 

/"* dt' 

6 i =5 + c i (x)-Ad i (x) J ^— (181) 

3 3 

The Ci and di are the two isocurvature-type (5 = 0, thus <5$ = 0) solutions. It happens that the relatively decaying 
isocurvature-type solution, i.e., c^-mode, temporally behaves the same as the peculiar velocity in the background, i.e., 
Di-mode. The relatively growing isocurvature-type solution Cj does not contribute to the u^, see Eq. (|172j) . From 
Eqs. p79| . p80| we have 

m-u=-[Vdi(x)+D i (x)-D(x)], (183) 
a 

which simply decays; the and D solutions are divergence-free and decoupled from the density perturbation, and 
are the peculiar velocity perturbation present in the background world model. 



with 



C. Comoving gauge 



To the linear order, only the hypersurface condition v — (the comoving temporal gauge) allows the density 
perturbation equation presented in the Newtonian form. Thus, we take 



v = 0, 



even to the second order. We take 7 = in Eq. (|141[) as the spatial gauge condition. Equation (|146[) gives 



E 



V U) V U)<* 



A" 





0)1/3 



Using this, Eqs. (|Tg2 > ([T4"5 |) give 

6-k= ~5,«X a + S K + \hY, ^0)^)%> + 3H E ^7f A_1V « 





0)1/3 



\ Uh + SvGg + c 2 ^) E W(3) v (3> +Uh + c 



1 " 





0)1/3 



Hi) 



K + C 



1 



\a a 



(i + %)H) 



3 „^ „ „„V^ MO) 



/3 

0)1/3 



[a (1 + %)) » (i)a ] ■ = --^ (v W/3 X ,/3 ) iQ + «« Wa - ^ («(0a«fi)) 



I. '3 



c 

"2^ 



E^o) v « («&)« 



0)/3 



3 





0)1/3 



(184) 
(185) 

(186) 



(187) 



(188) 



(189) 
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From Eqs. ( (fTM)) Jl87|) . and Eqs. (fT87l) - p89| we can derive, respectively, 



a 2 (S + ^6, a x a )\ - inGgS (1 + 5) = -\n^x a + ~/s a 



1 / A 



3 V C a 2X 



1 ,A 



2H + AnGg + ) E WS)«U)a + ( 6 ^ + c2 5 ) E ^^Va («8)«f 



a) I/? 



(c \ 1 ' c 4 



1 / A 

3lV : 



(k + p) ,Q v (i)a +2(n-tp) vf i)lc 





0)1/3 



(190) 



(191) 



Notice the 0(v^sV^j) a ) correction terms are present in Eqs. (|186|1 . (|187|) . (|190p even in the single component situation. 
Except for these 0(v^v^ a ) terms, the remaining parts of these equations coincide with the ones in the single 
component situation. 



D. Newtonian correspondence 



For K — 0, to the linear order, we have 

A A 

Xv = X - av = -av x , k v = -c— x v = c— v x , v {i)v = - v = w (i)x - v x , ip v = 0. 



(192) 



To the linear order we identify 



u = — cVv v , Ui = —cVv, 



Now, to the second order, we attempt identifying the Newtonian perturbation variables S, Si, u, and as 

«„ = — V-u, Xv = -u, ueVu, vu) v = - (iij - u) , 5 = S V , Si=Sa) v . 
a c w c 



(193) 



(194) 
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Using these identifications Eqs. (| 186|) - (| 1 9 1 1) can be written as 



(195) 



5+-\7-u= --V- (Su) + Hj2—\\i\^ -u| 2 + 3A- x V- [(u, ■ - u) V- (u 3 ■ - u 

Qj (X Q C I Zi 

3 

iv • (u + ffu) + 4ttG^ = ~V • (u • Vu) - C (t)Q/3 Qu |/s + 

+ Ef^{^ ( 4 ^" c2 ^) l u i- u l 2 + ( 127rG ^-c 2 ^) A^v-iK-tOv-K-u)]}, (196) 



j< + • u 4 - -iy • (SiUi) + \ [tyV ■ (m - u) - (u< - u) ■ Vy> + 2 « - u") 1 ^ C<g 



a a a 



(197) 



(198) 



+^i| Ui - u| 2 + 3if £ 77 ~ u > 2 + A ~' v ' K u i - u ) v ' ( u ^ - u )]| 

j 

-V • (in + Hm) + A-kGqS = -iv • (vu ■ Vu, ; ) - (~u a]f3 + C% 

Qj CL \ Qj 

+^GJ2 Qi^> \jVh - u l 2 + 3A ~ ly ' [("J - u ) V • K - u )]} > 
^ (aH) ' - AnGgd = ~ [aV • (<5u)]' + ^ V • (u • Vu) + ( H Ua|/J + C% 

- E { ( 47rG ^ - 77) l u J - u l 2 + ( 24 ^ - ^7) A ~' V • K u i - u ) V • K - u)]} , (199) 
-1 (a 2 ^)' - = -1 [aV • (^m)]" + \v ■ (u 4 • Vu,) + C® (-uf 3 + C«<* 

Q \ / Q, Qj \ Qj 



+ \ {A [u ■ (ui - u)] - V • [(iii - u) • Vu + u • V (ui - u)]} 



a' 
4irGg 



u, - u| 2 - SttG^T {\ Uj ~ u| 2 + 3A- X V • [(uj - u) V • ( Uj - u)]} . (200) 



c 



In Sec. I VII B II we have shown that, to the linear order, (u, — u) simply decays (oc a ) in an expanding phase. 
From Eq. (| 183|) we have 

- u = - [Vdi(x) + D<(x) - D(x)] , 

^Qjdj-Q, V • D = = V • T>i. (201) 

j 

Thus, Ui — u simply decays in an expanding background. If we ignore these contributions from the velocity differences, 
except for the presence of the tensor-type perturbation, Eq. (|200p coincides exactly with the zero-pressure limit of 
Newtonian result in Eq. (fH))) . Terms in the last two lines of Eq. (|200[) are relativistic correction terms which vanish 
for a single component case leading to Eq. (|199p : the second line in Eq. (|199[) also vanishes in the single component 
case. 

Ignoring quadratic combination of (u^ — u) terms, we have 

5 +-V -u = --\7 ■ (5u) , (202) 
a a 

iv • (u + Hu) + AttGqS = ~V-(u.Vu)- C^ af3 + C^j , (203) 

^ (a 2 6) ' - 4nGgS = ~ [oV • (Su)\ + \ V • (u • Vu) + (-u a \ + cj$) , (204) 

Qj V J Q Q V J 
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and 



1_ 1_ ,„ s 1 



a 



Si + -V • u, = --V • (Sim) + - 2pV • (Ui - u) - (u* - u) • Vv» + 2 « - u Q ) l/5 C<*2 , (205) 



a/3 



-V • (u, + Hm) + 4itGq5 = ~V ■ (m ■ Vm) - C ( ' )q/3 f -u a]f3 + C^ ) , (206) 
1 (a 2 ^) ' - 4vrG^ - ~ [aV • (^u,)]' + • (u, • Vm) + C*« Q< l/3 + C^ a A 

+ -- {A [u • (m - u)] - V ■ [(m - u) • Vu + u • V (u, - u)]} . (207) 

Equations (|202p - (|204p coincide with the density and velocity perturbation equations of a single component medium 
[1 21 ] ; thus, except for the contribution from gravitational waves, these equations coincide with ones in the Newtonian 
context. 

If we further ignore (iij — u) terms appearing in the pure second-order combinations, Eqs. (I205p - (|207[) become 

6 i + -V'U i = --V'(5 i u i ), (208) 
a a 

-V ■ (in + Hm) + AnGgS = -i V • (m • Vu,) - C^ af3 (~u al/3 + C%l) , (209) 
a a z \a 'J 

i (a 2 <5 4 ) ' - A-kGqS = ~ [oV • + • (u, • Vu.) + C*« Q< 1/3 + . (210) 



Notice that, by ignoring i-indices, Eqs. (|208|) - (|210[) coincide with Eqs. (|202|) - (|204p . In this context, except for the 
contribution from gravitational waves, the above equations coincide exactly with ones in the Newtonian context even 
in the multi-component case; compare with Eqs. (f8|).(|9"j). (|16p without pressure. In the single component situation such 
a relativistic/Newtonian correspondence to the second order was shown in [lj, [l2| . In the present case, the same 
equation valid in the single component is now valid in the multi-component case for the collective fluid variables. This 
justifies our identifications of Newtonian perturbation variables in Eq. (|194p . 

VIII. EFFECTS OF CURVATURE 

We consider a single zero-pressure, irrotational fluid. We take the temporal comoving gauge (v = 0) and the spatial 
7 = gauge. In the presence of background curvature the basic equations are presented in Eqs. (| 1 1 5[) - (| 1 1 7|) for 
nonvanishing pressure, or Eqs. Q188[) - (|191[) for zero-pressure multiple component fluids. By setting pressures equal to 
zero in Eqs. (JTT5J) - (jTTTJ) , or from Eqs. (fl87| - (fT9Q|) . we have 

S-k = --^S a r a +Sn, (211) 
or 

k + 2Hn- AttGqS = ~K, aX ' a + l -n 2 - ~ (c^ X ) + {c {t)a0 + ^X' q|/3 ) {c% + JrX,^) , (212) 

5 + 2HS - AitGqS = AttGqS 2 - 4 (<W*)' ~ ^,aX ,a + ^ 2 ~ I (A* 

+ (C^ + ± X ^) (C% + ± XMP ) , (213) 
A + 3K 



a* 

X + K = 0, (214) 

where Eq. (|214p is valid to the linear order. Compared with the situation with vanishing curvature, the effects of 
curvature in the above perturbed set of equations appear only in the linear-order relation between k and x m Eq. 

A. Newtonian correspondence 



Considering the successful Newtonian correspondence to the linear order even in the presence of the background 
curvature, we assume the identification in Eq. (|120p is valid to the second order. Then, to the linear order, from Eq. 
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(gTH) we have 



where u = Vu and x — Xv 



-V-u 



A + 3A 
-c - — x, 



a a a 

-av x . Thus, to the linear order, formally we have 

3A 



-x 



1 



u. 



(215) 



(216) 



a \ A + 3A, 

In the presence of curvature, the scalar-type perturbation can be handled by solving Eqs. (|21ip . (l212l) . (|214l) together 
with the identifications made above. We can formally separate the effects of pure curvature contribution. Using Eqs. 
([2I3,([216]), Eqs. (l2TT1).(l2T2l).(l2T?li become 
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From Eq. (|74p . the if term can be written as 
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IX. EFFECTS OF VECTOR-TYPE PERTURBATION 



The spatial C-gauge sets 

7 = = CM. (220) 

The remaining variables under this gauge condition are completely free of the spatial gauge modes and have unique 
spatially gauge-invariant counterparts. If we simultaneously take any temporal gauge which also removes the temporal 
gauge mode completely, all the remaining variables have corresponding unique gauge-invariant counterparts. The 
above statements are true to all orders in perturbations, see Sec. VI of From Eq. ([77)) we have 

0=- X , B«HfW- (221) 
a 

Thus, we have 

B a = - a X, a + ^\ C a p = vg§+C®. (222) 
As the temporal comoving gauge we set 

v = 0. (223) 

As we mentioned, the remaining variables under these gauge conditions are completely free of the gauge modes and 
have unique gauge-invariant counterparts to all orders in perturbation, see [Til ] . 



A. Linear perturbations 



To the linear order, the three types of perturbations decouple, and evolve independently. The rotational perturbation 
is described by Eqs. (206)-(209) in [ll|: 



A + 2A T M 8irG . \ (v) n 

■*L ) + — (/*+P) V a =0, 
C 4 



2a 2 
a 4 ( fi + p) v, 



(v) 



A + 2K n&° 



a 4 (jjl + p) 



2a 2 n+p' 



(224) 

(225) 
(226) 
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where the last equation follows from the first two. Compared with Bardeen's notation in [lj|, we have 

= c^-iff^QM ^) = *q«^(s( 1 )-^4 1 ))q«, 

«W = « C Q« - (>> - B«) Q« t£> + *W = <W - " ffi) Q« (227) 
where Bardeen's and Ws 1 ' , thus our u£ and v^ + ^a\ 

are related to the vorticity and the shear, respectively. 

From Eq. (|60p . to the linear order, we have 



0) 

uj al3 =av [aW , 



cr a p = a 



c. 



(t)/ 



(228) 



Bardeen called a 'frame-dragging potential'. The difference between and vi is crucially important to show 
Mach's principle including the linear order rotational perturbation in [27j ■ 

In the absence of anisotropic stress which can act as a sink or source of the angular momentum, we have 



Angular momentum oc a 4 (/i + p) u„ oc a 2 $^ oc constant in time. 
Thus, for vanishing anisotropic stress, we have 



oc 4 



In the zero-pressure limit we have 



a a z 



(229) 
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B. Second-order perturbations 



We consider a single component situation with general pressure. We set K = 0. To the linear-order we use 

A 

X = Xv = ~av x , k = k v =c—v x . 

a 

The scalar- type perturbation is described by Eqs. ([57]) , (jUI))) which give 

6 + 3H {c 2 s - w) 6 + 3H- -(1+w)k 
A 4 
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« + 2i?K - (fy + 3<5p) = - \3H + c 2 4 
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(232) 



(233) 



+ (Sn + 36p) a + ±Ha 2 + (a + 2<p) c 2 ^a + %a> a {a-<p) a + %C^a a m + (jjl + p) v^ a v^ 
cr 2 a z a z > a z c z 

+ (V* - (c 2 ^ x ) ^ + - ^ + (C*J - + . (234) 
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The vector-type perturbation is described by Eq. (|68|) which gives 



a a 5fi 

a ' 



(/i + p) «w] + - (/i + p) a, a + - (s P , a + nf ,g) 
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= p,cA a . 

We have L4 Q ] = From this we have 



2 A. 



(/i + p) a + 5p + -— n = a//A _1 V • A, 
3 or 



a 4 



A q -(A- x V-A) q 



Equation for the tensor-type perturbation (gravitational waves) follows from Eqs. (f6"6")) , (|102[) . 



(235) 

(236) 
(237) 



C. Zero-pressure case 



In the zero-pressure limit, Eqs. (|236|) . (|235[) give 
a = aA _1 V • A, 



aA n — — 
2 



(v) 



thus, a is purely second-order, and 
Equations (I233)) . ([234]) . (f237|) give 
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D. Newtonian correspondence 



In order to compare with Newtonian equations we continue identifying k as in Eq. (|120[) to the second order. To 
the linear order we identify 



u 



(244) 
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thus 



>) 



We introduce the following notations 



(245) 



(246) 



As mentioned below Eq. I|227j) . to the linear order, w a and C/ a are related to the vorticity and the shear, respectively. 
Equations PU J) - ([235 )1 become 



-V • u 
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V ■ (<5U) 
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247) 



(248) 
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From Eq. (|249|) we notice that the tensor-type perturbation does not affect the vector-type perturbation to the second 
order. The pure scalar-type perturbation also cannot generate the vector-type perturbation to the second-order; the 
same is true in the Newtonian case, see below Eq. (jX5j . 
From Eq. (|224[) . to the linear order, we have 



» 



1 / ck 



6\aH 



(250) 



where A: is a comoving wavenumber with A = — fc 2 , thus [fc] = 1. Since (ck) / (aH) ~ (visual- horizon)/ (scale), we have 



far inside horizon : 
far outside horizon 



U : 
U 



u = Vit - 

U = U a 



U ~ Vu, 



(251) 



Apparently, contributions of vector-type perturbation to the second order depend on the visual-horizon scale. 



Far inside the horizon, we can ignore c^a compared with In the matter dominated era we have 



5 5 a 2 AitGq ....... 



aH 

kc 



(252) 



where we used Eqs. (293), (329), (330) of [TJ. Thus, the third term in the right-hand-side of Eq. (|2^7| is (a-ff/fcc) 2 -order 
smaller than the first term. The fourth (and last) term in the right-hand-side of Eq. (|247[) is (aiy/fcc)[uW/(c#)]-order 
smaller than the first term. The third (and last) term in the right-hand-side of Eq. (|248p is also (aH /kc) 2 -order 
smaller than the first term. Thus, Eqs. 1(237 ]) -([2"39" ]1 give 
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a a a |p 

u + Hu) + AitGq5 = -^rV • (u • Vu) - C {t)ap 
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Hu^ = — [u • Vu - VA _1 V • (u • Vu)] . 
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(253) 
(254) 
(255) 



Thus, if we could ignore the tensor-type combination in Eqs. (I253[) . (|254[) . Eqs. (j253f) - (|255[) coincide exactly with the 
Newtonian equations: see Eqs. . (| 1 2 [) . (pT5|) ignoring the pressure terms and the subindices i. The vector-tensor 
combinations in Eq. (|253[) . (|254p are new relativistic contributions of the vector- type perturbations; compare these two 
equations with Eqs. (jTTJ) , (jT5]) which are valid in the absence of the vector-type perturbations. Notice the form of last 



term u 



MP 



c a\f) m fl^54[) which subtly differs from the expression u a \p in Eq. (TT8 



Contributions from the vector-type perturbation become more complicated near and outside the horizon scale. 

(v) 

The presence of vector-type metric perturbation Ufa , the scalar-type curvature perturbation tp, and the tensor-type 
perturbation C^i coupled with the vector-type perturbation give additional effects. 
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E. Pure vector-type perturbations 



In Sec. VII-E of [Til ] we have considered a situation with pure vector-type perturbation. As the analysis was made 
based on the fluid quantities in the normal frame, in the following we present the case based on the fluid quantities 
in the energy frame. Considering only the vector-type perturbation of a fluid, Eq. (|68p gives 



(v) 



c(A + 2K)IL£ 
2a 2 (pi + p) a 



Thus, for 11^ = we have 
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a^fi+p) 



« 4 (j* + pM 
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This differs from Eq. (365) of [llT ] which is due to the difference in the frame choice. The momentum constraint 
equation in Eq. (101) of [llj becomes 



2a 2 

Under the gauge transformation, from Eq. (|61[) we have 



c 4 a^HP) 1 



(258) 



(259) 



To the linear order, from Eq. (230) of [ll| we have Ba^ = B^ + 1;^' ■ We consider a gauge transformation from the 

C-gauge (Ca = 0j without hat) to the B-gauge (B^ = 0, with hat). We have Ba = 0, and B^ |c-gau g c 
Thus, 



(v) _ 



B (v) 



a I C— gauge 



*7 = -a a *W [% 



where we used B. 



a | C— gauge 



Thus, Eq. (|2351) gives 



u a IS -gauge 
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X. EQUATIONS WITH FIELDS 



A. A minimally coupled scalar field 



Equations in the case of a minimally coupled scalar field are presented in Eqs. (112)-(114) of [llj. The equation 
of motion in Eq. (112) and the full Einstein's equations in Eqs. (99)-(105) expressed using the normal-frame fluid 
quantities together with the normal- frame fluid quantities for the scalar field in Eqs. (114) all in [ll[ provide a 
complete set of equations we need to the second-order. The fluid quantities in Eq. (114) of [llj are presented in the 
normal-frame four-vector and it is convenient to know the conventionally used fluid quantities which are based on the 
energy-frame four-vector. These latter quantities can be read from Eqs. (88), (114) of [ll[ and Eq. (fT9"|) as 



Q N W = _i 
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+ P W ) (V<P -B a + AB a + 2V^C aP 
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5pW = S P N ^ 
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3<?' 



(262) 
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Thus, using Eq. (114) of [11| we have 



a<j) 
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nS = o. 



(263) 



Notice that no anisotropic stress is caused by a minimally coupled scalar field even to the second order in perturbations. 
The uniform-field gauge takes 8(f> = as a temporal gauge (slicing) condition to the second-order in perturbation. 
The uniform-field gauge gives = which is the comoving gauge, and vice versa. Thus, 







We also have 



• V W = o. 
Sn w - 8p w = 2V^8cb + V M &<?, 



and under the uniform-field gauge we have 



A(l-2A) + ^B a B a 



Equation (|263[) apparently shows that the vector-type perturbation Va' V ^ does not vanish to the second-order. 
However, the second-order quantities in right-hand-side depend on the temporal gauge condition for the scalar-type 
perturbations, and trivially vanish for the uniform-field gauge. We can also show that it vanishes for the uniform- 
density gauge or the uniform-pressure gauge where we have 5<f) — 4>A oc 8<f> to the linear order. In fact, we can show 
that the vector-type perturbation is not sourced by the scalar field to the second order by evaluating the rotational 
tensor uj a /3 in Eq. (|60p for the scalar field: i.e., for a minimally coupled scalar field we have 



(264) 
(265) 

(266) 



0. 



(267) 



to the second order. Thus, we conclude that a minimally coupled scalar field does not contribute to the rotational 
perturbation to the second order in perturbations. In fact, we can show that a single scalar field do not support 
vector-type perturbations to all orders in perturbations. As we take the energy frame, thus q a = 0, from Eq. (23) of 
fill we can show 



(268) 



Using the definition of the vorticity tensor in Eq. (I33|) we can show that u> a f, = 0. 

Using the fluid quantities in Eq. (|263p we can handle the scalar field using our non-ideal fluid formulation. The 
fluid equations in the energy frame, like Eqs. (|62p - (l68[) . remain valid in the case of scalar field with the fluid quantities 
expressed as in Eq. (|263|1 . The anisotropic stress vanishes and the entropic perturbation e is given as 



3 W = SpM 
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V, 



-(<!>) 



(<t>) 



<t> + v 4 



Under the comoving gauge v = 0, using Eq. (I266j) we have 

= 8^\ e(« = (l-^)^), 



(269) 



(270) 



which is a well known relation, now valid to the second order in perturbations. A fluid formulation of the scalar field 
to the linear order is presented in [UH^. Using Eq. (|270|) together with vanishing anisotropic stress, Eqs. (I62ll - (f68|) 
or (f88 |) - ([94|) provide the fluid formulation for a minimally coupled scalar field to the second order in perturbations. 
The perturbed equation of motion of the scalar field is presented in Eq. (112) of [TT[ | . 
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B. Minimally coupled scalar fields 



In the case of multiple minimally coupled scalar fields, the equation of motions in Eq. (119) and the full Einstein's 
equations in Eqs. (99)-(105) expressed using the normal-frame fluid quantities together with the normal-frame fluid 
quantities for the scalar field in Eqs. (121) all in [ll| provide a complete set of equations we need to the second-order. 
The fluid quantities in Eq. (121) of [ll| are presented in the normal-frame four-vector and it is convenient to know 
the conventionally used fluid quantities which are based on the energy-frame four-vector. These latter quantities can 
be read from Eqs. (121) of [ll[ and Eq. fl58|) as 



.m - 




(n) k 



N(4>) _ 



JV(0) 



8pM = 5 P 



ll af3 — L1 a0 



(271) 

where the normal-frame fluid quantities are presented in Eq. (121) of [ill] . Thus, by moving into the energy-frame 
fluid quantities from the normal-frame ones we have rather complicated terms which do not cancel out nicely any term 



in the normal-frame quantities. In single field case we had such cancelations, for example we have IT 



a.0 



in Eq. 



(|263|) . But, we do not have such a luxury in the multi-component situation. Apparently Tl^g in the energy- frame does 
not vanish and looks more complicated. Thus, in the multi-field situation we had better use both the fluid quantities 
and Einstein's equations all expressed in the normal frame: these are Eqs. (99)-(105),(121) in [111 ]. 

As the temporal gauge condition we can set any one field perturbation, say the specific l-t\i one 5<f>rf\, equal to zero 
which might be called the uniform-^) gauge to the second order. This apparently differs from the comoving gauge 

which sets = 0. In the multi-component situation we cannot take a gauge condition which makes v^'^ = 0. 
Thus, the multiple scalar fields source the vector-type perturbation to the second order in perturbations. The scalar-, 
vector-, and tensor- type decomposition of the anisotropic stress nj^g can be read by using decomposition formulae in 
Eq. (177) of [H|. 

In multiple-field situation, it is ad hoc and cumbersome (if not impossible) to introduce individual fluid quantity 
for each field variable even to the background and linear order perturbations, see [30j . 



C. Generalized gravity case 

In Sec. IV. D of [ll|we presented the equation of motion and effective fluid quantities in a class of generalized gravity 
theories together with additional presence of fluids and fields to the second order. The equation of motion is in Eq. 
(128) of [11|, and the full Einstein's equations in Eqs. (99)-(105) expressed using the normal-frame fluid quantities 
together with the normal- frame effective fluid quantities in Eq. (130) of [l l f p rovide a complete set of equations 
we need to the second-order. The effective fluid quantities in Eq. (130) of [111 ] are presented in the normal-frame 
four-vector and using Eq. (88) in [ll| we can easily derive the effective fluid quantities based on the energy-frame 
four-vector. As in the multiple field case in a previous subsection, by moving into the energy-frame, the effective fluid 
quantities become more complicated compared with the ones in the normal-frame. Thus, in these class of generalized 
gravity theories we had better use both the fluid quantities and Einstein's equations all expressed in the normal frame: 
these are Eqs. (99)-(105),(130) in [ll[. 
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XI. CURVATURE PERTURBATIONS AND LARGE-SCALE CONSERVATIONS 



In the large-scale limit the spatial curvature perturbation ip in several different gauge conditions is known to remain 
constant in expanding phase. Often the conservation properties are shown based on the first time derivative of the 
curvature perturbation. In order to show the conservation properties properly we have to construct the closed form 
second-order differential equations for the curvature perturbation. In the following we will derive such first-order and 
second-order differential equations for (p v , tp x , (fi K , and tp$. First we will derive equations for the linear perturbation 
including the background curvature and non-ideal fluid properties. Then we will derive equations for the second 
order perturbation assuming a flat background; including the background curvature is trivial, though. We consider a 
single-component fluid. 



A. Linear-order equations 



We introduce a combination 
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A-kG + p) 



This combination was first introduced by Field and Shepley in |28| |. From Eqs. 
Eqs. (HE]),®, and Eqs. flHEJ , dES]) , JH]) , respectively, we can derive 
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These equations were presented in |24j, 131] 

We can derive closed form second-order differential equations for Lp x , ip$, and (p K 
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90]) , dHj) , (HU) , 
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(274) 
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(276) 
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Equations ([277)) . ([278)) follow by combining Eqs. ([273)) . ([274)) . Equations ([279)) and ([280)) follow from Eqs. 
(EH)),([89l),([9l]),([92l), and Eqs. ® , (ES|l , (SB , © , respectively. 
From Eqs. l[59")) . ll9"0")) we can derive 
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These relations were presented in [171 ]. In the large-scale limit and in near flat background, thus ignoring A and K 
terms, we have 



$ ~ <p v ~ (p s ~ ip K , 



(282) 
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to the leading order in the large-scale expansion. 

In the large-scale limit, ignoring A terms, in near flat background and for an ideal fluid case, thus setting K = 
and e = = II, Eqs. (pTT jl -pSQ ]) give 

H^x) a H 2 ■ ( 283 ) 

tps <x ip K <x (285) 

Notice that if we simply ignore the A terms in Eqs. (I274j) - (|276j) we simply have ip v = ip$ = ip K = 0. In such a way 
we cannot recover the terms in the right-hand-side of Eqs. (|284j) . (|285p : these terms lead to decaying solutions (in an 
expanding era) in the large-scale limit and are higher order in the large-scale expansion compared with the decaying 
solution of ip x , see below. From Eqs. (|283j) - (|285|) we have general large-scale asymptotic solutions 

Vx = 4ttGC(x)- f a ^ p) dt + d(x)-, (286) 
a J H z a 

^CM + ^dwfjj^dt, (287) 

<p s = <p K = C(x) + |d(x) J J, (288) 

where C(x) and d(x) are integration constants which correspond to the relatively growing and decaying modes, 
respectively, in an expanding phase; in a collapsing phase the roles are reversed. The coefficients are fixed using 
the relations in Eqs. (|273p . (|274p . (|28ip . Notice that for the C-mode the relation in Eq. (|282p is satisfied, and simply 
remain constant. For the d-mode, tp v , ips, and ip K are A/ (aif ) 2 -order higher compared with the d-mode of ip x . The 
is one of the well known conserved quantity in the large-scale even in the context of generalized gravity theories 

In order to evaluate the solutions to the second order in the next section, we need complete sets of linear order 
solutions for different gauge conditions. For an ideal fluid, and for a minimally coupled scalar field such complete sets 
of solutions are presented in tabular forms in [2iJ l3l| . In the following we summarize such sets of solutions in an ideal 
fluid case for four different gauge conditions. From Table 8 of [3l| we have 

(JJ ft \ 2 jj ft 1 ft 

1 / adt I , S Y = = — 2C / adt, v y = —C—^ \ adt; 

a J J 3 if aH J x a 2 J 

^ = C, H Xv = C^j adt, ^ = _iip atJ = -?-A_C7(l-f J adt^, | = --A_cf J adt , 

H tig A 1 A H ( H 

(f6 = C, H Xs = C— adt, —=3a s = jTH C ^ v s = o^TEl C ?— 77 I 1 / adt 

a J H a z H z 6a z H z aH \ a J 

<p K = C, H Xk = C- ( adt, S H = -3^^a H = ~-^C, v K = ~-^-C^ [ adt. (289) 
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l + 3c 2 3a 2 iT 2 3 a 2 H 



For corresponding sets of solutions for a minimally coupled scalar field, see Table 1 of [29|. Compared with the notation 
used in [31| we have 7^ = —8TTG(/j,+p)av. The lower bounds of integration of solutions in Eq. (|289[) give behaviors of 
d-modes. For solutions without integration, the d-modes are A/(a7i) 2 -order higher than the non- vanishing d-mode, 
for example, see the solutions in Table 2 of [3l|. Thus, the d-modes are 

, H x 2k x ,1 
? x = -* x = d-, 6 x = --- = 2d—, v x =d-; 
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— , H X s = d, T7 = "q^T72 d / T' as = a^nr2 d "TTv ~ / ~3 ' Vs = q^77 d; 

a J an 6a z H z J a 6 9a z H z \a 3 H J a- 3 J 3 a 4 H 

A M ^ [* dt tt H ar o 1 + w 2 A 2 /•* dt 1 A 
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B. Second-order equations 



We assume K = 0. From Eqs. ([88 p . ([90 ]) . t|92 |> . Eqs. (JHH|) - d90J) , ((941) , Eqs. (88]) , (J93J) , and Eqs. (88]) , (89]) , (J9TJ) , respec- 
tively, we can derive 
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(293) 
(294) 



The perturbed order variables in Eqs. (|293p , (|294p are evaluated in the uniform-density gauge (S = 0), and the 
uniform-expansion gauge (n = 0), respectively. Equation (|293[) also follows from Eq. (4*1) evaluated to the second 
order. 

We can derive closed form second-order differential equations for tp v , ip x , (p$, and ip K 
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Equations (2^51) . (2^t5|) follow by combining Eqs. (pMj) . (j"2"9""2"j) . Equation (j2"9"7]) follows from Eqs. 
Equation (2"9"5|) follows from Eqs. (55]) , & , (H]) , (Ml) ■ The perturbed order variables in Eqs. (2SZ]) , (MS) are evaluated 
in the uniform-density gauge (5 = 0), and the uniform-expansion gauge (k = 0), respectively. 
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C. Large-scale solutions 

Now, we assume an ideal fluid, thus set e = = II. In the large-scale limit, thus ignoring the A/(aff) 2 -order higher 
terms, Eqs. (|256 > (|25g | give 

" I ( "° - m) + aHUe + ivGif+p) H U1 + Hn2 ) K 4«Gfr+pW (2 " } 

*P8 + -i- ( -rh x + ^Hm - Hn 3 ) - \n Q oc (300) 

/l + 3c 2 \ 1 1 

^ " 3^ l^^" 1 + "V " 3 n ° 00 (301) 

where the perturbed order variables in Eq. (|299[) are evaluated in the comoving gauge (v = 0). We already used the 
behavior of linear order solutions in Eqs. (j289[) . (|290[) in order to show that the right-hand-side of Eqs. (|296p - (|298[) 
vanish. Using the solutions in Eqs. (|289[) . (|290p we can show that 

(<p v -<£)' + 0(AC 2 , A 2 ri 2 ) cx ( ^ 2 J ffl3 , (302) 
(tps -<p 2 s)' + 0(AC 2 ,A 2 d 2 ) oc (<p K -ipiy + 0(AC 2 ,A 2 d 2 ) oc -1. (303) 
Thus, we have general large-scale asymptotic solutions 



^-rt=CM + ±dtof jj^dt, (304) 



.2 



^ - ^ = 9k - tp£ = C(x) + -d(x) y ^, (305) 

where C(x) and ci(x) are integration constants now including the second-order contributions, i.e., C = + C^ 2 \ 
etc. Ignoring the transient solutions in an expanding phase we have 

ip v = tps = ip K = C(x), (306) 

even to the second order in perturbations in the large-scale limit. 



XII. DISCUSSION 



In this work we presented pure general relativistic effects of second-order perturbations in Fricdmann cosmological 
world model. In our previous work we have shown that to the second-order perturbations, the density and velocity 
perturbation equations of general relativistic zero-pressure, irrotational, single-component fluid in a flat background 
coincide exactly with the ones known in Newton's theory, [12j . We also have shown the effect of gravitational waves to 
the second-order, and pure general relativistic correction terms appearing in the third-order perturbations, [l2l. [i"3j. 
Here, we presented results of second-order perturbations relaxing all the assumptions made in our previous work 
in [12j . We derived the general relativistic correction terms arising due to (i) pressure, (ii) multi-component, (hi) 
background curvature, and (iv) rotation. Wc also presented a general proof of large-scale conserved behaviors of 
curvature perturbation variable in several gauge conditions, now to the second order. 

Effects of pressure can be found in Eqs. (|123[) - (|125|) . As we emphasized, the effect of pressure is generically 
relativistic even in the background world model and the linear order perturbations. Still, our equations show the 
pure general relativistic effects of pressure (including stresses) appearing in the second-order perturbations. Effects 
of multi-component fluids can be found in Eqs. (|195|) - (|200|) . Although these equations apparently show deviations 
from Newtonian situation, in Sec. lVIlTJl we showed that if we ignore purely decaying terms in an expanding phase the 
equations are effectively the same as in the Newtonian situation. Effects of background spatial curvature K can be 
read from Eqs. l[2T7) l . ([215) 1 or Eqs. ()2TT )) - ([211) 1 . Effects of vector-type perturbation can be read from Eqs. ([247) 1 -([235 ]) . 
In the small-scale limit we showed that, if we ignore the tensor- type perturbation, the equations coincide with the 
Newtonian ones. 

Our results may have important practical implications in cosmology and the large-scale structure formation. Our 
new result showing relativistic/Newtonian correspondence in the zero-pressure irrotational multi-component fluids 
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is practically relevant in currently favored cosmology where baryon and dark matter are two important ingredients 
of the current matter content in addition to the cosmological constant. All equations in our work are valid in the 
presence of the cosmological constant. A related important result is the relativistic/Newtonian correspondence valid 
in the presence of rotational perturbation far inside horizon. Thus, inside the horizon scale, even in the presence of 
rotational perturbations we can still rely on the Newtonian equations to handle quasi-linear evolution of large-scale 
structures. As the spatial curvature in the present cosmological era is known to be small [33|, the possible presence of 
small spatial curvature may not be important in the second-order perturbations. Still, while the Newtonian equation 
is exactly valid to the linear order even in the presence of the spatial curvature, we have nontrivial general relativistic 
correction terms present to the second order in perturbations. Our second-order perturbation equations in the presence 
of pressure may have an interesting role as we approach early stage of universe where the effect of radiation becomes 
important. The importance of pressure to the second-order perturbations, of course, depends on whether nonlinear 
effects are significant in the early evolution stage of the large-scale structure during the radiation era and in the early 
matter dominated era. Realistic estimations of the diverse pure general relativistic contributions using the complete 
set of equations presented in this work are left for future investigations. 

In an accompanying paper we will investigate the effects of third-order perturbations of zero-pressure irrotational 
multi-component fluids in a flat background. This is one obvious remaining issue in our series of investigation of 
nonlinear cosmological perturbations where nontrivial general relativistic effects are expected. In the case of a single 
fluid we presented the pure general relativistic effects appearing in the third order in [131 ]. Corresponding results in 
the case of multi-component will be presented in [34j . 
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